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We study thermal and a transverse magnetic field response of a vortex line array confined to a 
plane with a single columnar pinning defect. By integrating out "bulk" degrees of freedom away 
from the columnar defect we reduce this two-dimensional problem to a one-dimensional one, localized 
on the defect and exhibiting a long-range elasticity along the defect. We show that as a function 
of temperature, for a magnetic field aligned with the defect this system exhibits a one-dimensional 
analog of a roughening transition, with a low-temperature "smooth" phase corresponding to a vortex 
array pinned by the defect, and a high-temperature "rough" phase in which at long scales thermal 
fluctuations effectively average away pinning by the defect. We also find that in the low-temperature 
pinned phase, the vortex lattice tilt response to a transverse magnetic field proceeds via a soliton 
proliferation "transition", governed by an integrable sine-Hilbert equation and analogous to the 
well-known commensurate-incommensurate transition in sine-Gordon systems. The distinguishing 
feature here is the long-range nature of the one- dimensional elasticity, leading to a logarithmic 
soliton energy and interaction. We predict the transverse-field — temperature phase diagram and 
discuss extension of our results to a bulk vortex array in the presence of a dilute concentration of 
columnar defects. 



I. INTRODUCTION 
A. Background and motivation 

The discovery of high-temperature superconductors al- 
most 20 years ago, among other things, has rekindled in- 
terest in the magnetic field [H] - temperature (T) phase 
diagram of type II superconductors^, generating a vigor- 
ous scientific activity. As a result, much has been clari- 
fied about the nature of equilibrium and nonequilibrium 
properties of vortex states in the presence of thermal 
fluctuations, pinning disorder and electrical ("super"-) 
current^^^, leading to a rich phase diagram. In particu- 
lar, in contrast to a mean-field phase diagram, thermal 
fluctuations drive a first-order melting of a vortex lattice 
over a large portion of the phase diagram into a resistive 
(although with large conductivity and diamagnetic re- 
sponse) vortex liquid^^^'', that, from the symmetry point 
of view is qualitatively identical to the normal state. Fur- 
thermore, as was first shown by Larkin^^d2^ arbitrarily 
weak pinning disorder, on sufficiently long Larkin scale 
(that diverges in the limit of vanishing disorder) always 
disrupts translational order of the vortex lattice. It has 
been argued theoretically^' with a limited experimen- 
tal support^'*, that in the resulting state, vortices are 
collectively pinned into a vortex glass characterized by 
an Edwards- Anderson^^ like order parameter, exhibiting 
a vanishing linear mobility and therefore a vanishing lin- 
ear resistivity. 

While the original proposal for the vortex glass state 
was made in a context of intrinsic, short-range corre- 
lated (point) disorder, it was soon appreciated that the 
beneficiary effects of pinning can be enhanced by intro- 
ducing artificial pinning centers by, for example, elec- 
tron and/or heavy ion irradiation, with the latter re- 



sulting in a forest of columnar pinning defects, that is a 
particularly effective pinning mechanism^^. The result- 
ing anisotropic vortex glass was dubbed a Bose-glass"'^^'^^ 
because of its mathematical connection with interact- 
ing two-dimensional (2D) quantum bosons pinned by 
a quenched (time- independent) random 2D potential 
This connection allowed understanding many of the prop- 
erties of the anisotropic vortex glass from the correspond- 
ing quantum Bose-glass phase. 

One key feature of the anisotropic vortex glass that dis- 
tinguishes it from the corresponding (putative) isotropic 
one is the existence of the "transverse" Meissner effect,^* 
namely a vanishing response to a field H± < 
applied transversely to columnar defects and vortex 
lines. This expulsion of the transverse flux density, 
B±, that has received considerable experimental^" and 
simulations^^ support, corresponds to an effectively di- 
vergent anisotropic vortex glass tilt modulus^^, that 
in the quantum correspondence maps onto a vanishing 
superfiuid density in the Bose-glass phase. The de- 
tailed theoretical description of the transverse Meiss- 
ner effect (as well as other properties of the phase) 
has been predominantly limited to noninteracting vor- 
tex lines""^®'^^. Although these are supported by scaling 
theories^^'^^'^^'^'' (borrowed from the variable-range hop- 
ping theory for electronic systems^^) that do incorporate 
effects of both disorder and interactions (clearly essen- 
tial for the very existence of the Bose-glass phase), with 
the exception of functional RG analysis^^, a detailed in- 
teracting description is limited to simulations. ^^'^^ This 
is not surprising, as a description of strongly interacting 
random systems is a notoriously difficult (with few excep- 
tions) unsolved problem, whose solution is at the heart of 
understanding many of the interesting condensed matter 
phenomena. 
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One way to incorporate strong interactions is to ap- 
proach the problem from the vortex solid (rather than 
the vortex liquid starting point) pinned by a random po- 
tential. Potential difficulties with this approach are a 
proper incorporation of topological defects (dislocations 
and disclinations) that tend to proliferate in the presence 
of quenched disorder and external perturbations. 

Recent analytical real-space renormalization group 
(RG) study^^ that demonstrated stability of a 3D weakly 
disordered random-field XY model to a proliferation of 
topological defects (vortices), provide a strong argument 
for the stability of an elastically disordered but topolog- 
ically ordered vortex Bragg-glass phase postulated and 
studied in detail by Giamarchi and Le DoussaP^. These 
studies therefore give support to treatments of vortex 
solids^^'^®"^^ that ignore the notoriously difficult-to-treat 
topological defects. Furthermore, even if Bragg-glass 
is unstable to dislocations, for weak disorder disloca- 
tions will be dilute, with physics on scales smaller than 
their spacing expected to be well-described by the vortex 
Bragg-glass phenomenology. ^ 3 , 28-3 1 

A suppression of topological defects can furthermore be 
facilitated by a planar confinement of vortices, realized 
in layered high-Tc superconductors (e.g., BISCCO)'^^ or 
by artificially prepared multilayers^^'^^, where for a mag- 
netic field directed along the planes, vortices are well lo- 
calized to 2D. The resulting planar (2D) vortex array 
pinned by point disorder, where dislocations are excluded 
by construction is in fact the "toy" model studied by 
Matthew Fisher^^, that motivated his original proposal 
of a vortex glass phase in bulk superconductors. 

Motivated by these ingredients, in this paper we study 
a (l-l-l)-dimensional vortex array confined to a planar 
slab of thickness w in a presence of a single planar colum- 
nar defect, illustrated in Fig.l. This system was intro- 
duced and first studied in great detail in Ref. 35. 




FIG. 1: A (l+l)-dimcnsional (planar) vortex lattice pinned 
by a single columnar defect studied in this paper. 

As should be clear from the above discussion, such 
a "toy" model'^^ may be relevant to the regime of far - 
separated (by d, compared to vortex spacing y/^o/B) 
columnar defects, accessible for flux density far exceed- 
ing the columnar-defect matching field = (j)o / cP 



where 0o = hc/2e « 2.1 x 10"'^ G -cm^ is a fundamental 
quantum of flux. As first investigated in Ref. 35, we study 
the response of such a (l-l-l)-dimensional vortex array to 
a planar tilting magnetic field H±, applied transversely 
to the columnar defect, as illustrated in Fig. 2. 
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FIG. 2: A (l-l-l)-dimensional (planar) vortex array induced 
by a magnetic field H applied in the plane of the slab, at an 
angle to a single columnar pinning defect. 

The rest of the paper is organized as follows. We con- 
clude the Introduction with a summary of our main re- 
sults and predictions, hcuristic;ally extended to a dilute 
concentration of columnar defects. In Sec. II we derive 
the appropriate (l-l-l)-dimensional continuum model for 
a single defect and discuss its ingredients. By integrating 
out "bulk" (away from the defect) degrees of freedom we 
reduce this model to a (O-l-l)-dimensional model confined 
to a defect, and characterized by a long-range elasticity 
along the defect. In Sec. Ill we study the effect of trans- 
verse magnetic field H±^ and demonstrate that tilting of 
the vortex lattice away from a columnar defect proceeds 
via a novel "commensurate-incommensurate transition" 
controlled by a proliferation of solitons. In Sec. IV we 
study effects of thermal fluctuations and demonstrate 
that this system exhibits a ID "roughening" -like tran- 
sition. We explore its consequences for the vortex po- 
sitional correlations, and construct a H±-T phase dia- 
gram. We conclude in Sec.V with an extension of these 
results to an experimentally relevant case of a dilute con- 
centration of columnar defects (allowing for genuine tran- 
sitions) and close in Sec.VI.with a summary of our study. 

B. Summciry of results 

The body of the paper is primarily devoted to the 
study a planar (l-l-l)-dimensional vortex array at a ID 
vortex density no = 1/a, in the presence of thermal fluc- 
tuations and induced by an external planar magnetic field 
H = HzZ + Hx applied at an angle to a single columnar 
defect'^^. As we will show below, even this "toy" prob- 
lem is quite rich, providing insight into the bulk (2-|-l)- 
dimensional multi-defect problem. It has the added ben- 
efit that it can be analyzed in detail analytically. This is 
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not surprising as such a planar classical vortex array is 
a cousin of a one-dimensional quantum problem, a Lut- 
tinger liquid, that is known to be exactly solvable and to 
exhibit rich phenomenology.^^ In fact our classical anal- 
ysis of the vortex problem has strong formal connections 
to the work of Kane and Fisher^^ who studied a Lut- 
tinger liquid in a presence of a single localized impurity, 
a problem that admits exact analysis^^. This connection 
was first emphasized and fruitfully utilized by Hoffstetter 
et al.^^ and Polkovnikov, et aL**", although we will not 
take advantage of it. As is well known^^, the Luttinger 
formalism is equivalent to the classical theory of vortex 
lattice elasticity that we (and Ref. 35) employ here. 

Our work has a strong overlap with that of Ref. 35, 
particularly on the finite temperature analysis for a van- 
ishing transverse field in an infinite single-pin system, 
as well as a large transverse field, where tilt response 
is analytic. Where this overlap exists our predictions 
are in complete agreement with those found in Ref. 35. 
However, our emphasis is on the low-temperature, strong 
coupling regime, where tilt response is highly nonlinear, 
and can only be understood in detail in terms of vortex 
lattice solitons, that proliferate at a novel commensurate- 
incommensurate crossover, that, we argue, turns into a 
genuine sharp phase transition for a dilute concentration 
of columnar defects. 

Hence, as we describe in more detail below, in the 
presence of a dilute concentration of columnar defects 
the planar vortex array exhibits two phases in the trans- 
verse field H± - temperature T phase diagram: a low 
T, H± "commensurate" pinned/ aligned (C) phase and 
a high T, H± "incommensurate" depinned/ tilted (I) 
phase. As illustrated in Fig. 3, a novel finite temper- 
ature commensurate-incommensurate phase transition 
separates the two phases. 

For a vanishing transverse field, the phase transition 
is akin to a thermal roughening transition^^^^"*, that is 
closely related to a zero-temperature superfluid-insulator 
transition in a resistively-shunted Josephson junction^^, 
the opaquc-to-transparcnt impurity transition in a Lut- 
tinger liquid"^*, and spin-boson and dissipative impurity 
models^^-^^. 

Before we discuss our results, a disclaimer is in or- 
der here. Because a single columnar defect cannot pos- 
sibly compete with bulk degrees of freedom, with its ef- 
fects vanishing in thermodynamic limit, we expect a two- 
dimensional free energy that is analytic, and therefore no 
true phase transition can take place in a 2D thermody- 
namic limit. However, as in, for example, the Kondo 
problem'*'', where one considers the effect of an impurity 
on the bulk electron gas, bulk effects are predicted only 
once a finite density of defects is considered. Here too 
the transitions that we discuss are for the elastic degrees 
of freedom localized on the columnar defect, that is a 
boundary critical phenomenon. Although the ID rough- 
ening transition is indeed a genuine one, the contribu- 
tion of the associated (nonanalytic) free energy to the 
bulk two-dimensional system vanishes in the 2D ther- 
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FIG. 3: H_i_ - T phase diagram illustrating the commensurate 
(C) and incommensurate (I) phases separated by a continu- 
ous CI phase transition at H]^{T) (full curve), where vortex 
lines in the vicinity of far-separated columnar defects tilt via 
proliferation of solitons. As the soliton density grows with 
increasing H± (dashed curve), the solitons overlap near the 
upper-critical transverse field H'^{T), beyond which the sys- 
tem crosses over to a smoothly tilted vortex lattice. The 
lower-critical field, H'^{T) vanishes with increasing pin sepa- 
ration d, eliminating the CI transition (but not the roughen- 
ning transition at Tp) for a single columnar defect. 



modynamic limit. Furthermore, the CI (vortex line tilt- 
ing) "transition" takes place at a lower-critical transverse 
field H'^{L), that, in the case of one pin is driven to 
zero in the thermodynamic limit (L oo). It only be- 
comes a genuine sharp phase transition for a finite density 
1/d of columnar defects, with the lower-critical trans- 
verse field H'j^{d) set by the columnar-pin spacing d.^® 
Since this latter case is the one most easily accessible ex- 
perimentally, we heuristically extend our rigorous single 
pin results to a finite dilute concentration of indepen- 
dent columnar defects. This being a notoriously difficult 
unsolved problem, we expect this heuristic extension to 
break down on sufficiently long scales, where collective 
pinning effects become important. ^'''"^^'^^"^^ 

Our residts can be divided into two categories: thermal 
effects at a vanishing transverse field, i.e., for the mag- 
netic field (and therefore induced planar vortex array) 
aligned with columnar defects, and a low-temperature 
response of the vortex array to a tilting (transverse) mag- 
netic field. 

For a magnetic field aligned with columnar defc^cts we 
show that a vortex array undergoes a "roughening" -like 
transition at 



(1.1) 



between a low-temperature state in which each colum- 
nar defect effectively pins the vortex lattice around the 
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pin for T < Tp and a phase in which its pinning effects, 
even near a columnar defect, vanish at long scales. In 
above, K and B are vortex lattice tilt and compressional 
moduli, whose vortex density and dependence on other 
parameters (e.g., vortex interaction) can be computed 
from microscopics^^'^'^''*^'^^. This roughening-like transi- 
tion is related to a one-dimensional long-range interact- 
ing — jp) Ising model**^, that is well-known from 
a number of other physical contexts, most prominently 
the Kondo problem**^ , and more recently, in a dissipative 
Josephson junction^^ and a Luttinger liquid in the pres- 
ence of an impurity^^ . In relation to the latter work, we 
note that unlike the quantum case where there is no sim- 
ple way to tune the Luttinger parameter^^, here it can 
be simply tuned by temperature and vortex density.^^ 

As is usually the case for such (e.g., roughening''^"^'^, 
Kosterlitz-Thoulcss^^, and other topological) transitions 
there is no local order parameter and phases are distin- 
guished by long-scale behavior of correlation functions. 
The low-temperature T < Tp, pinned phase is character- 
ized by a finite mean-squared vortex fluctuations at the 
location of the defect (that we take to be a; = 0), with 
correlations away from the defect given by^^ 



{u{z, x)u{0, x')) 



K{\x\ + \x'\f + Bz-' 



K{x - x'f -h Bz"^ 



(1.2) 



On scales longer the pinning length ^ (defined in 
Eq.4.10, below), the corresponding average density ex- 
hibits Friedel-like oscillations^^ given by 



{n{x,z))o 



no + 2nG ( 1 cos(27rnoa;), (1.3) 



where a = Max [a, ^^B/K] and 

27r ksT 



7] 



(1.4) 



In contrast, at high temperature T > Tp, vortex ther- 
mal fluctuations effectively average away the effects of 
the pin, leading to phonon correlations that diverge log- 
arithmically with sample size. The connected phonon 
correlation function is finite and is given by 



{{u{x,z)-u{0,0)f) 



In 



B 



(1.5) 



Despite the irrelevance (in the RG sense of the term) 
of the columnar defect for T > Tp, the average density 
also displays perturbative Friedel oscillations^^ 



ncav fW f a \ ^ ^ ,^ ■, 
i(x,z)) « "-0 + -^^Y B I Id ' cos(27rnoa;), 



(1.6) 



with a stronger power- law exponent than that for T < Tp 
and an amplitude that vanishes with the strength of the 
pinning potential v. 

The low-temperature T < Tp, the pinned (commen- 
surate) phase is distinguished from the rough (incom- 
mensurate) phase by a transverse Meissner response to 
a magnetic field H± applied transversely to the colum- 
nar pin. Namely, we find that for a field smaller than a 



lower-critical transverse field 



W 



iln^ 

In 



' K 2TTd 



(1.7) 



where L is the length of the sample along the columnar 
defect {z), w is the slab thickness, and d is the columnar 
pin spacing. At low temperature ^ « is approximately 
given by 



^0 



27r2 



For < vortex lines in the wide vicinity 



1/2 



(1.8) 



(1.9) 



of a columnar pin remain aligned with it, therefore ex- 
hibiting a bulk transverse Meissner effect for a < A" 
spaced array of pins. Of course, (as for the thermally- 
driven depinning transition discussed above, here too), 
because a single pin cannot compete with the bulk mag- 
netic energy, away from the defect beyond this screening 
length the vortex lattice is always aligned along the 
applied magnetic field. Related to this, for a single pin 
(c? ^ oo) the critical transverse field ff^, Eq.1.7 clearly 
vanishes in the thermodynamic L oo limit. 

For a transverse field stronger than a continuous 
vortex lattice tilting transition takes place into a tilted 
(incommensurate) state. It proceeds via a proliferation 
of solitons, with soliton density ns{H±) and the average 
vortex tilt at the location of the columnar defect growing 
continuously beyond 



H < H^, 
H > H^^ , 



(1.10) 
(1.11) 



In the incommensurate phase, the screening length out 
to which the vortex lines are aligned along and pinned 
by the columnar defect diminishes with increasing soliton 
density and is given by 



A/. = 



1/2 



1 



B 

KJ 2Tms{Hs_y 
1 



(1.12) 
(1.13) 



As we show below, for finite transverse field H±^ , Friedel 
oscillations decay exponentially away from the columnar 
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defect with length also given by A^. We estimate the 
upper-critical transverse field H^, for which soliton lat- 
tice becomes dense, to be given by 



H. 



c2 



(1.14) 



In the presence of fluctuations, Hc2 field marks a 
crossover from a nonlinear soliton tilted regime to a uni- 
formly tilted state. For large H± exceeding H^, the 
vortex lattice tilts smoothly and the screening length re- 
duces to A/i ~ 1 /H-^ , in this limit coinciding with Friedel 
oscillations decay length found by Affleck, et al^^. 

The finite temperature CI phase boundary in Fig. 3 is 
given by 

H^,{T) « H^^{0)^1-^^ (1.15) 

The corresponding upper-critical crossover field H^{T) 
(dashed curve in Fig. 3) for weak pinning (^o ^ o,\jK/B) 
can be estimated by using Hc\{T) and ^(T) inside 
Eq.1.14, 




T/\Tj,-T\ 

(1..16) 



For strong pinning (^o a^jKjB), _ff^(T) is given by 
the above expression, but with set equal to a^jK/B. 

In the remainder of the paper we demonstrate results 
summarized above. 



II. MODEL 

A. Vortex lattice elasticity with a transverse 
magnetic field 

In a type-II superconductor, for fields above a lower- 
critical field magnetic flux penetrates in a form of inter- 
acting vortex flux tubes, with average density determined 
by the applied magnetic field. ^ At low temperature and in 
the absence of disorder, a periodic array (Abrikosov lat- 
tice) of repulsive, elastic vortex lines forms, whose elastic 
description'*^ can be derived from the Ginzburg-Landau 
theory for the superconducting order parameter, that it- 
self, under certain conditions, is derivable from the micro- 
scopic theory of superconductivity. As usual, transcend- 
ing such a detailed derivation, on sufficiently long length 
scale the elastic vortex lattice energy functional can be 
deduced purely on symmetry grounds. In d-dimensions, 
it is formulated in terms of a — 1-dimensional Eulerian 
phonon field (Goldstonc mode of the spontaneously bro- 
ken translational symmetry) u(x, 2) = x — Xj describing 
a transverse vortex lattice distortion at a d-dimensional 
position (x, z) relative to a perfect vortex array charac- 
terized by (xj, z). 



A planar vortex array, that we take to be confined 
to the X — z plane, is characterized by a scalar phonon 
field u{x, z), describing x-directed vortex distortion with 
a continuum elastic Hamiltonian given by 

We; = ^ y dxdz [K{d^u - hf + B{d^uf] , (2.1) 

where K and B are tilt and compressional elastic mod- 
uli, respectively, that we take to be phcnomcnological 
parameters. ^° The parameter h encodes the effect of an 
additional magnetic field H±, applied transversely to 
the columnar defect (z-) axis, with h = H ^/ {(p^^n^) = 
Hj_/Hz in a sample that (other than the pin) we take for 
simplicity to be isotropic. 

Vortex pinning, characterized by a weak potential 
Vpin{x, z) can be easily incorporated through its coupling 
to the local vortex density ny{x, z) via 



(2.2) 



Tip = I dxdzVpin{x, z)ny{x, z). 



As with any periodic elastic medium, vortex density is 
given by 

ny{x,z)^no-nod,u + J2'^G, e*°''(^+"(^'^», (2.3) 



with vortex lattice distortion u(x, z) entering through the 
variation of the long-scale density fluctuation, —riodxU, 
and via the variation of the phase, GpU of the vor- 
tex density wave given by the last term. In above, 
Gp = 27rnop = 27rp/a (a the vortex lattice constant, 
p G Z) spans a one-dimensional reciprocal (to x) lat- 
tice and no = 1/a the average x-projected vortex ID 
density. Above representation for ny{x,z) can be de- 
rived in a standard way from its microscopic definition 
ny{x,z) = ^iS{x — Xi{z)) in terms of vortex-line con- 
figurations Xi{z), by the use of the Poisson summation 
formula'^^, with the key periodic (last) term arising from 
vortex discreteness. 

For the problem of a single ^-directed columnar defect, 
we can approximate the pinning potential by an attrac- 
tive zero-range form V^in = ~VocS{x), with Vq and c its 
effective strength and range, respectively. For simplicity, 
and without loss of qualitative generality, we include only 
the lowest harmonic, characterized by the minimal recip- 
rocal lattice constants G±i = ±G = ±27r/a to model the 
periodic vortex density. Furthermore, by minimizing the 
total energy, it is easy to show that the long-range part 
of the density variation, —nodxU has a simple effect of a 
small shift in vortex positions, Uo{x, z) = — ^"^'^'^ sgn{x), 
that is constant and positive to the left of the defect 
and constant and negative to the right of the defect. It 
thereby slightly increases the average vortex density, but 
only at the location of the columnar defect, x = 0, that 
can be absorbed into the background density. Dropping 
an unimportant constant, the pinning Hamiltonian then 
reduces to 



dz cos(Gu(0, z)), 



(2.4) 
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localized at the defect at x = 0, with v = 2nGV(jC. The 
resulting total Hamiltonian, Ti, = Hei + 'Hp 

W = 1 y dxdz [K{d^u - hf + B{d^uf] 

-V j dz cos {Gu{Q,z)), (2.5) 

is reminiscent of the well-known sinc-Gordon model de- 
scribing a broad spectrum of commensurability phenom- 
ena in condensed matter physics, ranging from crys- 
tal surface roughness to topological defects in ordered 

media^^"''^ 

There is, however, an essential difference in that the 
nonlinear pinning term is localized at x = 0. As a re- 
sult, away from the defect, the system is harmonic and 
therefore solvable by elementary methods. Clearly, as il- 
lustrated in Fig. 2 away from the defect the vortex lattice 
must asymptote to that of a columnar defect-free config- 
uration, that simply follows the transverse field 



^{x, z) = hz, for X — > ±oo. 



(2.6) 



obtained by minimizing H for v = 0. As we will show 
below, we can fruitfully take advantage of the locality of 
the pinning potential by "integrating out" (eliminating) 
the bulk elastic degrees of freedom away from the defect, 
thereby reducing the two-dimensional problem to an ef- 
fective one-dimensional nonlinear one, that can be solved 
exactly. The approach is quite similar to Kane-Fisher's 
treatment of a point impurity in a one- dimensional elec- 
tron liquid.^* 



B. Reduction to one-dimensional model 

As noted above, because of the short-range nature of 
the pinning potential in 7i, Eq.2.5, vortex degrees of free- 
dom, u{x, z) away from the columnar defect at a; = 0, 
are governed by a harmonic Hamiltonian. As a result, 
they are simply related to vortex distortion at the colum- 
nar defect, allowing us to eliminate u{x, z) in favor of 
u(0, z). To automatically satisfy the boundary conditions 
dzu{x,z)\z=o,L = h, induced by finite h, it is convenient 
to shift to a new phonon variable, measuring vortex lat- 
tice distortion relative to Uoo, 

u{x,z) = u{x,z) — hz, (2.7) 

in terms of which the Hamiltonian becomes 

n = dxdz [K{d,uf + B{d^uf] 

-V j dzcos[G{u{0,z) + hz)]. (2.8) 

At zero temperature a reduction to an effectively one- 
dimensional model can be most straightforwardly done 
by solving the Euler-Lagrange equation for u{x, z) for a 



prescribed (but arbitrary) distortion u(0, z) = U{){z) on 
the columnar defect (In Appendix A, we present a com- 
plementary constrained functional integral-based deriva- 
tion of above result, that extends to finite temperature) 

{Kdl + Bdl)u{x, z) = Buo{z)dJ{x). (2.9) 

The local stress term (source term on the right-hand-side) 
is chosen so as to produce a vortex lattice distortion at 
a; = to automatically satisfy the boundary condition 
u{0,z) = uo{z). Standard Fourier analysis leads to a 
solution 

ii(x,5,)=uo(g.)e-(^/^)'''l'''ll^l, (2.10) 

which, when substituted into the Hamiltonian, H, Eq.2.8, 
and integrating over x reduces to 

Ho = j ^|g,||fio(g,)|2 

-V j dzcos[G{uo{z) + hz)]. (2.11) 

The two-dimensional nature of the underlying vortex lat- 
tice is captured by the nonanalytic form {\qz\) of the 
effective one-dimensional elasticity in Hq. As in other 
examples of a low-dimensional system coupled to a bulk 
system (e.g., a crack or a crystal surface in a bulk solid^*) 
it encodes long-range interactions of one-dimensional de- 
formations mediated through bulk (away from the colum- 
nar defect) degrees of freedom, as can be easily seen by 
reexpressing Ho in terms of uo{z) 



Ho 



27r 



■// 



dzdz' 



I ( uo{z) - uo{z') 



z — z' 



-V j dzcos[G{uo{z) + hz)]. (2.12) 

The long-range elasticity qualitatively distinguishes 
this system from a one-dimensional sine-Gordon model 
characterized by short-range and therefore analytic {q^) 
long-scale elasticity. As we will see in Sec. IV, the asso- 
ciated enhanced stiffening of elastic distortions is what 
allows this one-dimensional system to undergo a finite 
temperature roughening phase transition, in contrast to 
a one-dimensional sine-Gordon model. For reasons that 
will become clear below and by the analogy with the sine- 
Gordon model, we refer to this system as the sine-Hilbert 
model. In the next section we will study the sine-Hilbert 
model at zero temperature but finite tilting field h, in or- 
der to characterize a vortex lattice response to a magnetic 
field H±_ applied transversely to the columnar defect. 



III. ZERO-TEMPERATURE TRANSVERSE 
FIELD RESPONSE 

A. Bulk tilt response 

For /i = much is known about the sine-Hilbert 
model, Hq, Eqs.2.11,2.12, as it arises in many different 
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physical contexts including rcsistivcly-shuntcd Josephson 
junctions"'^ and Luttinger liquid transport in the presence 
of an impurity^^. In the present context, these findings 
relate to a vortex lattice at finite temperature, a study 
that we will undertake in Sec. IV. Furthermore, using in- 
verse scattering transform Santini, Ablowitz and Fokas^^ 
have shown that the classical sine-Hilbert model is inte- 
grable and admits soliton solutions, a finding that we will 
make use of below. 

Less is known about the finite h phenomenology, the 
study of which is facilitated by returning to the Uo{z) = 
Uo{z) + hz displacement field. In terms of ^0(2;), the 
Hamiltonian becomes: 



Ho = 



-V J dzcos[Guo]- 



uo{z) - uo{z') - h{z - z') 



(3.1) 



For a finite transverse field, a good starting reference 
point are two competing solutions 



uc{z) = 0, 
ui{z) = hz, 



(3.2) 
(3.3) 



that minimize Holuoiz)] in the /i — > and h ^ 00 
limits, respectively. By analogy with commensurate- 
incommensurate (CI) phase transitions*^'^'^ we refer to 
the corresponding phases as "commensurate" and "in- 
commensurate" , respectively. In the commensurate 
state, uc{z), vortex lattice aligns with the columnar de- 
fect, minimizing the pinning energy, while raising the dia- 
magnetic energy. In the incommensurate state, Ui{z), 
instead, the vortex lattice aligns with the external field, 
thereby ignoring the defect and sacrificing its attractive 
pinning energy. For a sample of extent L along z, the 
corresponding total energies in the two cases are given 
by 



Ec = 
Et = 0. 



h^L'' - vL, 



2lT 



(3.4) 
(3.5) 



The extensive scaling, ^ of the magnetic energy in 
Ec is expected from Eq.2.10, that shows that elastic dis- 
tortion on scale L along z, imposed at a; = (in the com- 
mensurate phase corresponding to misalignment with ap- 
plied field) decays over length L{B/ KY^^ into the bulk, 
leading to a region of area ~ L-^{B/Ky^-^ with a finite 
diamagnetic energy cost. As discussed in the Introduc- 
tion, in the thermodynamic limit this bulk diamagnetic 
energy always dominates over the linear-with L pinning 
energy, and vortex system is in the incommensurate state 
for arbitrarily small transverse field h. However, for a fi- 
nite L, we find from Eqs.3.4, 3.5, that for a single pin, the 
tilting transition between uc and uj states takes place 
at the critical ("thermodynamic bulk") tilting field 



he 



27TV 1 



(3.6) 



that, as expected, vanishes in the thermodynamic limit. 
L — !■ CX3. For a finite density of columnar defects spaced 
by d, it is clear that the bulk critical field saturates at a 
finite d-dependent value given by 



> 1/2 

27TO ^ ' 1 

Li/2 ■ 



VBK 



< d 



L\/§ > d 



(3.7) 



This oversimplified picture of the tilting transition, ig- 
noring solitons (see below) is summarized in Fig. 4. 







E(h) 



h 



Commensurate / Incommensurate 



-vi 



FIG. 4: Oversimplified (ignoring solitons) expression for the 
energy E{h) of a vortex lattice in the presence of a transverse 
field h, indicating the aligned (C) to tilted (I) transition at 
he- 



B. Tilt-solitons 

Two states, uc, uj, are only appropriate in the h ^ 0, 
h ^ 00 limits, respectively. However, in analogy with 
other systems, where there is competition between elas- 
tic and pinning energies'*^"'*'', we expect and find (see be- 
low), that the tilting transition into an incommensurate 
state is driven by soliton proliferation transition above 
a lower-critical field hd, that preempts the bulk tran- 
sition at he found above. The soliton state above hci 
then continuously approaches the fully incommensurate 
ui{z) solution in the /i — » 00 limit, when solitons become 
dense. The existence of a lower-energy soliton solution 
can be seen by a simple inspection of the Hamiltonian 
Ho, Eq.3.1. It stems from the periodicity of the pinning 
energy, Hpin[uo] = Hpin[uo+a], that microscopically cor- 
responds to its independence of which of the identical 
vortex lines in the array is pinned by the columnar de- 
fect. A soliton at zo corresponds to a solution uosiz) that 
switches at zq between two adjacent minima of the peri- 
odic potential. As illustrated in Fig. 5, from the 2D per- 
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spcctivc such sohton describes a switching between two 
adjacent vortex Unes locahzed on the columnar defect. 



FIG. 5: A two-dimensional perspective of a ID soliton along z, 
corresponding of an exchange between two neighboring vortex 
lines localized on the columnar defect. 

At zero temperature the soliton shape is characterized 
by a single length scale, 



6 



V2^ 



(3.8) 



that can be read off from Ho by bala ncing the elastic 
and pinning energies, that scale like \/ KBc? and ?;^0) 
respectively. This soliton width is set by the balance 
between the pinning and elastic energies. The former 
(latter) is minimized by the most (least) abrupt shift be- 
tween the adjacent pinning minima, with reflecting 
this through its dependence on elastic moduli and pin- 
ning strength, increasing with B, K and decreasing with 
V. We note that the effective ID modulus is a geomet- 
rical mean of the vortex compressional (B) and tilt (K) 
moduli, as 2D distortion corresponding to the soliton con- 
figuration along the defect in Fig. 5 involves both tilt and 
compression of the 2D vortex lattice. 

It is convenient to express length scales along z in units 
of the soliton width, ^o, and trade in the displacement 
field uq[z) for a dimensionless phase field 



(2) = ^Mo(^;), 



(3.9) 



where throughout the paper we will used "hat'ed" sym- 
bols to denote dimensionless quantities. The Hamilto- 
nian then reduces to a dimensionless form 

2 



Wo 



1 

47r 



dzdz' 



(j){z) - <f>{z') - h{z - z') 



z — z' 



I 



dz cosli 



(3.10) 



where 

= ^ov, 

is the soliton energy scale and 



(3.11) 
(3.12) 

(3.13) 



is a dimensionless measure of a transverse field. 

These energy and length scales are sufficient to deter- 
mine all qualitative ingredients of the soliton- driven tilt- 
ing (CI) transition. However, integrability of the Euler- 
Lagrange equation 







{z)-4>{z') 



-l-sin^(i) =0, 
[z — z'Y 



(3.14) 



found by Santini, et al.^^ allows for a quantitatively exact 
analysis. In above, the integral must be defined as the 
principle value with the singular point z — z' excluded, 
and we have used free boundary conditions on (^(L) and 
0(0) in deriving Eq.3.14. These lead to a boundary con- 
dition 



5s < 



(3.15) 



that supplements Eq.3.14. It encodes the condition that 
at the edge of the sample vortex lines tilt to slope h to 
follow the external magnetic field. 

A single-soliton solution to the above sine-Hilbert 
equation, Eq.3.14 was discovered by Rudolf Pcicrls^" in 
his seminal study of an edge dislocation in a crystal, and 
later rediscovered and considerably extended (to multi- 
solitons, dynamics and proving integrability) using the 
inverse scattering method by Santini, et al.°^. It is illus- 
trated in Fig. 6 and given by 



,{z) = -2ArcTan- 



1 



zo 



(3.16) 



This solution can be verified by a direct substitution of 
<t>s(z) into Eq.3.14, using Hilbert transforms that we re- 
view in Appendix D. 

Substituting 0s(i) into Hq, Eq.3.10, and using Hilbert 
transforms (with details presented in Appendix E), we 
find the energy of a single soliton state 



= eoEc + eoEsi 



(3.17) 
(3.18) 



where dimensionless single soliton energy, Egi , computed 
in Appendix E is given by 

£^^1 = 27rln(^) -/i(2L-27r), (3.19) 

L = L/^o, and Ec the dimensionless energy of the com- 
mensurate state (cf. Eq.3.4) 



Er 



477 



h^L^ - L. 



(3.20) 
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FIG. 6: Soliton solution (j>s(z), Eq.3.16 to the sine-Hilbert 
equation, Eq.3.14. Inset shows dz4>s{z)- 

From Egi, Eq.3.19 it is clear that soliton energy be- 
comes negative for h > hd, with 

a 1 , el/ /„ X 

" (3.22) 

that, as expected (given that there is only a single pin) 
vanishes in the thermodynamic limit. As with other CI 

phase transitions'*^''*'', this leads to soliton proliferation 
at hci^^, that preempts the bulk transition at he, ap- 
proaching the fully tilted incommensurate state ui{z), 



only as h ^ hc2- The latter tilt field is defined by when 
the solitons begin to overlap, namely ng{hc2) ~ Co"^- 
However, in contrast to conventional CI transition, where 
the ratio of the lower critical field to the thermodynamic 
one IS of order 1 constant (= 4:/{V2tt) « 0.90 for the CI 
transition in the sine-Gordon model), here 

Vanishing of this ratio for large L demonstrates the im- 
portance of solitons in driving the vortex lattice tilting 
transitions. 



C. Tilt-soliton proliferation transition 

The vortex lattice tilting angle, 0{h), related to trans- 
verse flux density B±{H±), via b± = tan 9 = B^/B^ 
is determined by the soliton density ns{h) through the 
relation 

6_l(/i) = ans{K) (3.24) 

The soliton density ns{h) for h > hd in turn is deter- 
mined by the balance of the soliton chemical potential en- 
ergy h (that induces solitons) and the soliton repulsive in- 
teraction. In the dilute soliton limit (for hd < /i >C hc2), 
we can approximate the latter by a sum of pair-wise 
soliton interactions. This is determined by the energy 
£"2(^1, ^2) = eo-E'2(-Zi, -S2) of two one-solitons, localized at 
zi and Z2, separated by a large distance Z1 — Z2 '> 1, with 



-^2(^1, 22) = 'HQ[(f>sl + (t)s2], 



dzdz 



, ( (jisliz) + (j)s2{z) - (psl{z') - 'Ps2{z') - h{z - Z') 



z — z' 



(3.25) 

- j dzcos[(l)d{z) + (t>s2{z)]. (3.26) 



A straightforward calculation along the lines of the com- 
putation of El gives (with details presented in Appendix 
F): 

E2{zi,Z2) = Ec + '2Ed{h)+Vs{zi- Z2), (3.27) 
where the soliton two-body repulsive interaction is given 

by 



Vs{z) = 27rln 



52 -H 4 



4-77, 



(3.28) 



47rln— -F47r, for 1 < |5| «: L/2, (3.29) 

2\z\ 



and is illustrated in Fig. 7. 



The energy of the Ng soliton lattice is then given by 

(3.30) 



Ens = Ec + NsEd+J2^'^^^i~^i^' 

i<j 



Ec + 2L{hd - h)N, + -V,{L/2)N^ 



Ec + L^ 



2Chd-h)hs + -Vs{L/2)hl 



(3.31) 

, (3.32) 



where, because of the long-range nature (logarithmic) 
of the interaction Vs{z), the sum in the interaction en- 
ergy is dominated by the longest length scale, L, lead- 
ing to its iV, = n^L^ stronger-than-extensive growth. 
This is in qualitative contrast to standard CI T = 



10 



4Tiln(U4) 




FIG. 7: Soliton interaction potential Vs{z), illustrating a long- 
range logarithmic repulsion. 



transitions^^'^^'^^ , where the interaction is short-ranged 

and the corresponding sum is dominated by the small- 
est term ^5(0/^0); leading to an extensive interaction 
energy." 

Minimizing E^g[ng] over the soliton density fig = 
Ng/L, we find the advertised CI soHton proHferation 
transition^^ at hc\, with the soliton density 



hs{h) 



0' 

nso{h ■ 



hci), 



h < hci 
h > hci, 



(3.33) 



growing linearly with the transverse field h. hgo is 0(1) 
constant given by 



(3.35) 



figo = , (3.34) 

VgiL/2y ^ ' 

1 

The final numerical result, Eq.3.35, above is quite crude, 
providing only an order of magnitude estimate, as we 
have not carefully treated the case of a soliton near the 
edge of the system. 

The dimcnsionless soliton density fig (measured in 
units of 1/^0) increases to 0(1) when h = h^^^ exceeds 
hci by 1. This corresponds to the physical soliton density 
ns{hc2) ~ 1/Co of a dense lattice of overlapping solitons, 
with the upper-critical transverse field hc2 given by 

hc2 = hci + (3.36) 

consistent with physical picture illustrated in Fig. 16 and 

summarized by Fig. 8 

Using solution Ng{h) = hg{h)L, Eq.3.33 inside Ens 
we find the vortex lattice energy as a function of the 
transverse field h: 



in 
47r 



L, h < hci 

L- ^{h- hci)^P, hci<h^hc2. 

(3.37) 



Mp(h) 




Conmensurate 



Incommensurate 



FIG. 8: A sketch of soliton density, ns{h), and the number of 
vortex lines pinned per columnar defect, Mp{h), as functions 
of the transverse field h. 



As expected the negative soliton energy (last term) can- 
cels ofi^ the I/^ misalignment-field energy (first term), as 
the incommensurate state is approached with increasing 
soliton density in the large h ^ /td limit. "^"^ This energy 
and the corresponding phase diagram are illustrated in 
Fig.9. 







Coinmensurate Incommensur; 



E(h) 




-vL 



FIG. 9: Vortex lattice energy E{h) as a function of transverse 
field h, illustrating soliton-drivcn tilting CI transition at the 
lower-critical field hci preempting the bulk thermodynamic 
critical field he- At high upper-critical transverse field hc2 
the system crosses over to a fully discommensurated tilted 
state with a smooth tilt response. 
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D. Bulk vortex lattice distortions 

Having established the existence and computed the de- 
tails of the T = vortex array tilting transition in terms 
of the one-dimensional field uq{z), illustrated in Fig. 10, 
we now turn to the computation of the associated bulk 
vortex lattice distortion uo{x,z). The connection be- 



/ 



f 



FIG. 10: A sketch of the function uo(a) for an array of 5 
solitons. 




FIG. 11: A sketch of the function ut^(z), Eq.3.42, character- 
izing the periodic part of the soliton train solution wo- 



These respectively contribute to the oscillatory and uni- 
form part of the (negative) transverse magnetization. 

To compute the Fourier transform mqI'Zz); we extended 
/ieff-2 part of Mo(-z) beyond the system size L to a continu- 
ous periodic function with period 2L illustrated in Fig. 12. 




twccn ?io(^) and mo(.t,z) is established through a key 
relation, Eq.2.10 (here we use the subscript to denote 
T = field configurations, obtained by minimizing the 
energy) 



'"o(a;,gz) = 'ito(gz)e" 



-\AI>-h{iz) 



(3.38) 



that shows quite clearly that a one-dimensional distortion 

with a wavcvcctor along the defect penetrates into the 
bulk a; ^ over a screening length 



(3.39) 



A real-space one-dimensional distortion in an soli- 
ton state is given by 



^o(-2) = y^ns(-z - Zg) - hz, 



a=l 



uq{z) — (h — nsa)z, 



(3.40) 

(3.41) 



with Us{z) = {a/27r)(j)s{z) the single soliton solution from 
Eq.3.16. Physically dzUo{z) is proportional to transverse 
magnetization. In Eq.3.41 we defined a periodic part 
coming from the soliton array (illustrated in Fig. 11) 



"0(2) = ^Us{z- Za) - 



HsttZ, 



(3.42) 



a=l 



that oscillates around with period l/ris, and the re- 
maining misalignment part hesz, with an effective tilt- 
field 



hes = h-ns{h)a, 

= h — 2'7rnso\h — hci\ 



(3.43) 
(3.44) 



FIG. 12: A sketch of a tilted component of uo{z), periodically 
extended beyond the system size L. 



From this decomposition it is clear that the distortion 
Uo{z) appears on two characteristic wavelengths: the sys- 
tem size L along z (wc take the system to be infinite along 
x) and the soliton spacing nj^. Hence we predict that 
the associated bulk distortion, uq{x, z) will display an 
average misalignment with the applied field (i.e., on av- 
erage align with the columnar defect), that extends over 
a length 



AS = ^ 



1/2 



(3.45) 



along X and will exhibit soliton-induccd oscillations with 
wavelength ~ n^^ that penetrate over scale 



A. = 



KJ 2'Kns{h)' 



(3.46) 



into the bulk, x 7^ away from the columnar defect. 

Above qualitative discussion can be elevated to an ex- 
act calculation. Using Eq.3.38, the real-space bulk dis- 
tortion is given by 

uoix,z)=hz + J ^uo(3z)e-^^'*l«^ll^le^9»". (3.47) 



Putting together above ingredients, the Fourier trans- 
form un{qz) c;an be easily computed, and, when inserted 
into Eq.3.47 gives 
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(3.48) 



4 ■r-^ 

uo{x,z) = hz T^LheS y ^ 



i-iy 



- (2m +1)2 



1 _ (1 _ !^)e- 



i(2m+i)|a.| sin[^(2m + l)z] + -Y ie-\^''^"^fl^l sinpTrrisp^], 

(3.49) 



p=l 



TT 



Although no approximation is required to compute 
Mo(a;,2), Eq.3.47, in Eq.3.48 we have considerably sim- 
plified the Fourier transform by approximating the soli- 
ton part U{^(z) by a function — rigoz, periodically ex- 
tended with period nj^. This approximation is valid 
away from the dense soliton limit, i.e.. for // <C hci-, where 
riga <C 1. In going to Eq.3.49 we have furthermore sim- 
plified uq(x, z) by keeping only the first harmonic in each 
of the two terms and approximated the first sine by a line 
with a proper slope (determined by boundary conditions) 
valid ior z <^L away from the edges of the sample. The 
corresponding configurations of the vortex lattice in the 
commensurate state, h < hci, slightly inside the soli- 
ton state, h Ri /i^, intermediate regime < h <^ hc2, 
and deep in the soliton state h 3> hd are illustrated in 
Figs. 13, 14, 15, and 16 respectively. 

As anticipated by the qualitative discussion above, in- 
fiuencc of the columnar defect depends on the tilt field 
h and ranges from the maximum value A^, Eq.3.45 in 
the commensurate state, down to A^, Eq.3.46 in the soli- 
ton state. Upon increasing h beyond the lower-critical 
field hci, the soliton density increases and the penetra- 
tion length decreases according to 



Xh{h) = 



a fB_ 



1/2 



1 



\h - hc\\ 



(3.50) 



a 

2^ 



1/2 



\, for /ici <C /i <C hc2, 
h 





FIG. 13: Vortex lattice (T = 0) configuration uo{x,z) for a 
subcritical tilt field h < hd, showing expulsion of tilt over 
a region of size L X {B/Kf'^L around the pin, but with 
penetration of tilt on scale at x = growing with |a;| to 
full penetration beyond |,r| ~ {B/KY''L. 



(3.51) 




saturating at the microscopic lattice scale a as hc2 is 
approached. The size of the pinned vortex "cloud", 
Xh allows to define an important dimensionless number 

Mp(/l)35.64 



Mp(/l) = 



Xh 



(3.52) 



that gives the number of vortex lines effectively pinned 
by a single columnar defect. As illustrated in Fig. 8 it 
ranges from M™"*^ « {B/Ky/'^L/a {d/a for finite density 
of pins) in the commensurate (h < hci) state down to its 
minimum value of 1 for h hc2- 



FIG. 14: Vortex lattice (T = 0) configuration uo(x, z) for just- 
above critical tilt-field h = h^, with a single soliton along 

the pin. It shows the expulsion of tilt over a region L x 
{B/K f''^L/2 around the pin, with the tilt penetration length 
Xh2 ~ {B/ Ky^^ L/2 along \x\ reduced by a factor of 2 relative 
to the case of no solitons, illustrated in Fig. 13. 



IV. FINITE TEMPERATURE 
THERMODYNAMICS AND CORRELATIONS 

Having established the vortex lattice T = tilt- 
response to a transverse magnetic field, we now turn 

to the study of finite-temperature properties of the vor- 
tex lattice pinned by a single columnar defect. "^^ As 
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FIG. 15: Vortex lattice (T = 0) configuration uo{x,z) for 
tilt-field hci < ft <§C hc2, with throe solitons along the pin. It 
shows the expulsion of tilt over a region L x {B / Ky^^ L/4 
around the pin. with the tilt penetration length Xh4 = 
[B / KY^^ L/A along \x\ reduced by a factor of 4 relative to 
the case of no solitons, illustrated in Fig. 13. As discussed in 
the text, in general, the A'^-soliton penetration length is given 
by A;.jv = {^Y'''L/N, ^ Xn/Ns = n-'(B/Ky/^ 





















^^^^^^^ 



FIG. 16: Vortex lattice (T = 0) configuration ucix^z) for tilt- 
field h hc2, giving a dense solitori array and a vanishing 
tilt-expulsion length. A resulting vortex lattice tilt response is 
nearly that of a pin-free system, with the distortion confined 
to the immediate vicinity of the columnar defect. 

usual, the thermodynamics and corresponding correla- 
tion functions can be computed by integrating the vortex 
phonon configurations weighted by a Boltzmann factor, 
g-'H/ksT 1 2^ ^ with Z the corresponding partition function 
Z = /[du]e-«/'=«^. 



A. Thermal depinning transition for vanishing 
tilt-fleld, h = Q 

We first consider the case of a vanishing tilt-field, 

/i = 0, and show that the effective one-dimensional, long- 
range interacting Hamiltonian, T-Lq , to which our 2D vor- 
tex problem has been reduced, exhibits a thermal de- 
pinning transition, qualitatively similar to a well-known 
roughening transition in 2D systems. "^^^^^ 



The indication of the existence of the transition comes 
from a computation of the thermodynamics pcrturba- 
tively in powers of the pinning potential strength v. As 
for the 2D roughening (and related) transition, this per- 
turbation theory diverges (at long length scales) for low 
temperatures T < Tp even for an arbitrarily small v, 
but is convergent for high temperatures T > Tp, with 
the critical pinning temperature Tp defined below. As 
usual to make sense of the associated infrared diver- 
gences near and below Tp we employ the momentum-shell 
renormalization-group (RG) transformation^^ on the ID 
Hamiltonian Ho 

Ho = VkB J ^\q^\\Uo{q,)\^ - V J dzcos[Guo{z)]- 

(4.1) 

Namely, in the partition function for the model, we inte- 
grate out perturbatively in v short-scale phonon modes 

u^{qz) in a narrow shell GzC^^ < \qz\ < Gz aroimd 
the short scale cutoff Gz = G^B/K = {2TT/az), with 
^.^ <i; 1. This allows us to express the partition function 
in terms of integrals over the remaining effective course- 
grained phonon modes u^(gz), with \qz\ < GzC^'^^. In 
order to massage the resulting Hamiltonian into the orig- 
inal form Ho, i.e., to bring the reduced ultra-violet (UV) 
cutoff back up to Gz , we rescale lengths and wavevectors 
according to 

= q'ze-'', (4.2) 
z = z'e^^. (4.3) 

For convenience, to keep the vortex lattice reciprocal vec- 
tor G in the argument of the pinning potential fixed, 
we choose not to rescale the real-space phonon field 
Uq{z) = u'q{z'), which implies a rescaling 

u<{q)=e-''u',iq'), (4.4) 

for its Fourier transform. With these transformations the 
effective coarse-grained Hamiltonian returns to its origi- 
nal Ho form, but with effective, ^-dependent elastic and 
pinning parameters. As can be easily checked, to lowest 
order in v, the one-dimensional stiffness VBK remains 
unchanged, and the effective pinning strength trans- 
forms according to 

= ,,e(i-''/2)« ' (4.5) 

where t] is defined by root-mean-squared phonon fluctu- 
ations 

rj ^ ^{{u>{z)f)o, 

^ G^ ksT r^' dqz_ 
2TrS£ Jc^e-" \<lz\' 
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Treo 



(4.6) 



Because of the aforementioned relation of our model to 
a variety of other problems in condensed matter physics, 
this RG analysis is in fact quite familiar from those other 
contexts^^'^'*''*^'^^. As advertised it clearly predicts a 
phase transition at = 2,^^ corresponding to 



27reo, 



(4.7) 

(4.8) 



between the pinned phase for T < Tp and the dcpinncd 
phase for T > Tp. The two phases are distinguished 
by the relevance and irrelevance of the pinning potential 
V, respectively. More physically, for T > Tp thermal 
fluctuations of the vortex lattice at long scales effectively 
average away the effects of the pinning potential (but see 
below), reducing it relative to the elastic energy. 

In contrast for low temperatures T <Tp, as indicated 
by the RG flow, Eq.4.5, the strength of the pinning po- 
tential, no matter how weak at the lattice scale grows 
relative to and becomes stronger than the typical elastic 
energy. Quite clearly, since the computation of the RG 
flow in Eq.4.5 is done perturbatively in v the flow ter- 
minates on scale ^(T) = UzC^* when elastic and pinning 
energies become comparable. On scales longer than ^(T), 
the displacement field is well-localized at a minimum of 
the periodic pinning potential, which therefore can be 
safely Taylor-expanded to quadratic order and treated 
as a "mass" for uo{z). Hence ^ is determined by the bal- 
ance of the eff'ective pinning and elastic energies on scale 




-r,/2 



(4.9) 



with the effective pinning strength Vdt = ^(^/a^)^''/^ < 
V reduced by thermal fluctuations of modes in the range 
of scales Uz < z < ^. Solving Eq.4.9 for ^ we find 



^0, 



^0 > a\ 



(4.10) 



that, for weak short-scale pinning is, as expected, expo- 
nentially lengthened by thermal fluctuations 



B. Depinned phcise, T > Tp, h = 

Irrelevance of the pinning potential for T > Tp al- 
lows us to compute correlation functions in the dcpinncd 
phase perturbatively in v. To zeroth order, it is computed 
via a Gaussian integral with Hamiltonian Hei, Eq.2.1, at 
long scales giving 



Ct>t,{x,z) = {{u{x,z)-u{Q,0)f), 

,dqz l-e'^-- 



2kBT 
ksT . 



f dqxC 
J (2^ 



)2 Bql + Kqr 



+ ^z^ 



(4.12) 



(4.13) 



where we used vortex lattice spacing o as a natural short- 
scale cutoff on X. At the location of the defect at x = 
0, CT>Tp{0,z) can be equivalently computed using the 
effective ID Hamiltonian "Hq, Eq.2.11, to zeroth order 
{v = 0) giving 

Ct>tAO,z) = {{uo{z)-uo{0)f), (4.14) 

keT f^'dqz._^ , 
(1 - cosg^z). 



(4.15) 



Vortex density correlations can be similarly computed 
using Eq.2.3 and phonon correlations above. To zeroth 
order in v, the average density is given 

{n{x, z))o = no- no(a,u)o + 2nGRe(e'2-"o(x+«(x,.))^^^ 

■n/2 

^^0+ I t\/^ I 2nG cos(27rnoa;), 





no, 



(4.16) 



C(T)«6e'^^/(^''-^\ for 6> 




(4.11) 



with, a = \n{^o\/ B / K / a) > 0. 

We next turn to the computation of correlation func- 
tions in each of these phases. 



with fluctuations vanishing in the thermodynamic (i 
oo) limit. Physically we do not expect this to be the case, 
as the pinning potential breaks translational invariance in 
X, even when it is irrelevant. This reveals that for T > Tp 
the pinning potential is in fact dangerously irrelevant, as 
dropping it completely (as in above calculation) restores 
translational invariance and gives a result that is not even 
qualitatively correct for {n{x,z)). Indeed to correctly 
capture the behavior of (n(x, z)), we need to compute it 
to at least 1st order in v. The irrelevance of u for T > Tp 
guarantees the convergence of such perturbation theory. 
We find^s 

{n{x,z))-no « 2^Re / dz'(e^G(-+"(-'^)) cos[C?«(0, z')])o, 

kbT J 

^ cos Gx [ d^'e-*([«(-.^)-«(o.^')]^>o 
■n-i 



nova K ( a 
c„- — —\l-^[-ri] cos(27rnoa;). 



" ksT V B \ \x 



(4.17) 
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(cri = 0{1) dimensionless constant) showing that even 
in the phase where it is irrelevant, the pinning poten- 
tial leads to Friedel oscillations in the density, a result 
missed by simply setting v = 0. This illustrates often 
under-appreciated distinction between irrelevance in the 



RG sense and unimportance of an operator in the physi- 
cal sense. 

Density two-point correlations can also be straightfor- 
wardly computed. Using above phonon correlations and 
Eq.2.3, we find^^ 



/ / \ / f r\\\ I / W I f I r\\\ ^O^B^ I Bz (x X ^ c\ 1 ( ^ 1 /e\ \ 

{n{x,z)n{x ,Q)) « {n{x, z))o{n{x , 0))o + ^ ^ [^^(^ _ ^,^2 + g^2]2 + _ ^ b^2 j cos(27rnoar) 

(4.18) 

I 



with the first term given by Eq.4.17, and with the last 
term approximated by its lowest order (in v) translation- 
ally invariant expression. 

C. Pinned phase, T <Tp, h = 

In the low temperature pinned phase, v is relevant, 
growing with increasing length scale relative to the elas- 
tic energy. On scales longer than ^, we can therefore 
approximate its strong pinning effects by replacing it by 
a Dirichlet boundary condition on u{x, z) at the loca- 
tion of the pin away from its ends (for a more systematic 
justification of this, see Appendix B), namely taking 

w(0,z) = Mo(z) = 0. (4.19) 

Bulk phonon and density correlations with this boundary 
condition can still be computed exactly. As derived for a 
general case in Appendix B, a pinned phonon correlation 



function can be expressed in terms of a pin-free correla- 
tion function, i.e., from a purely harmonic elastic theory 
with w = 0: 

Gt<t^{x,x';z) = (u(x,z)u(.t',0)), (4.20) 
= Go{\x-x'\,z)-Go{\x\ + \x'\,z). 

(4.21) 

where Gq{x,z) is w = two-point phonon correlation 

bmction. As expected, because of the pin at x = 0, 
Gt<Tp(x, a;'; z) is clearly not translationally invariant in 
X, depending on both x and x' . It is easy to see that it 
vanishes identically for x and x' on opposite sides of the 
pin, e.g., a; > and x' < 0, showing that because of the 
pin such phonon fluctuations are completely uncorrelated 
on length scales longer than ^. On the other hand, for 
x and x' on the same side of the pin (e.g., a; > and 
x' > 0), above expression is given by 



Gt<Tp{x,x';z) = Go{\x — x'\, z) — Go{x + x' , z), for a; > 0, a;' > 0, 



-{[u{x, z) - u{-x', 0)]2)o - -{[u{x, z) - u{x', 0)]2)o, 



K{x + x'Y + Bz 



K{x - x'f -h Bz"^ 



(4.22) 
(4.23) 



r 



and at the same point [x — x' , z = z'), reduces to 

Gt<t,{x,x;0) = Go(0,0)-Go(2x,0), (4.24) 
= ^([n(a;,0)-«(-ar,0)Ho, (4.25) 



2ny/KB ^\ a 



(4.26) 



In above, we natur ally c ut off separation at short scales 
by a = Max[a,i^B/K], since the Dirichlet boundary 



condition, Eq.4.19 is only valid on scales longer than a. 
Utilizing this result, we can easily compute the average 
vortex density at point x, z 



{n{x, z)) « no + 2nG 



2\x\ 



cos(27rnox). (4.27) 



{n(x, z)) in the pinned phase also displays Friedel oscilla- 
tions similar to that of the depinned phase, Eq.4.16,4.17, 

but with an amplitude that is nonperturbativc in v and 
power-law fall-off with exponent 77/2, rather than 77 — 1. 
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The latter is continuous at the roughening transition, 
where rj{Tp)/2 = r]{Tp) — 1 = 1. This is in agreement 
with the result first found by Affleck, et al. by utilizing 
an equivalent Luttinger liquid phenomenology^^. 



27r In 



2L 



1 - 4 1 In 



2LV* Tp)"'^{T) 



(4.28) 



D. Finite T and finite tilt-field h 

The calculations of previous sections have focused on 
either thermal properties at a vanishing transverse field, 
or on finite transverse field response at a vanishing tem- 
perature. As wc show below, these can be extended to a 
general point T > 0, /i > on the phase diagram, Fig. 3. 



1. Finite T tilting transition 

The zero-temperature results for the tilting transition 
of Sec. Ill can be extended to a finite temperature by 
utilizing RG analysis and matching from Sec. IV A. As 
in a related sine-Gordon commensurate-incommensurate 
transitions^^'^^ the basic nature of the transition re- 
mains the same but with thermally renormalized effec- 
tive parameters. The extension of the zero-temperature 
transverse field boundary, Eq.3.22 to a finite-temperature 
phase boundary hci{T) is determined by the difference in 
free- energies of the commensurate and a single-soliton 
states, computed for T < Tp. The former is sim- 
ply given by the sum of the energy of the commensu- 
rate state Ec and the entropic free-energy contribution 
Fu = — fcsTlnO„ of the phonon modes. 

The computation of a single-soliton free energy is 
more complicated, since at a low temperature the pin- 
ning potential is relevant and therefore requires a strong- 
coupling analysis. This can be done by matching free- 
energy calculations on scales shorter and longer than the 
strong coupling correlation length ^ (T) . As is clear from 
the analysis leading to ^(T), Eq.4.10, for ^ ciz (corre- 
sponding to weak pinning v) the pinning potential is sub- 
dominant on scales < z < £,{T) and the free-energy is 
dominated by due to Gaussian phonon fluctuations. 
On scales longer than £,{T), the pinning potential domi- 
nates and free-energy is determined by the sum of three 
contributions: (i) a single soliton energy Ei{T), with 
short-scale cutoff ^(T) replacing in Eq.3.19, (ii) soli- 
ton's positional entropic contribution —kBT\n{L/^{T)), 
and (iii) free-energy contribution due to phonon fluc- 
tuations about a single soliton configuration. Noting 
that phonon fluctuations about the commensurate and 
single soliton states are approximately the same, i.e., 
Fu ~ F^ + F^, we find that the free energy difference 
between a single soliton and the commensurate state is 
given by: 



This leads to the prediction for hc\{T) quoted in the In- 
troduction: 



2L 



(4.29) 
(4.30) 



In going from the first to second line we used a 
temperature-independent approximation for ^(T), 
Eq.4.10, valid for large L. 



2. Correlations atT > 0, h> 

Correlation functions at a finite h and T can also be 
computed. The nature of valid approximation and their 
resulting form strongly depends on three possible regimes 
on the phase diagram, Fig. 3. 

(a) Incommensurate (fully-tilted) state, h > hc2{T): 

The strongly incommensurate, high tilt and tempera- 
ture regime h > hc2 (T) , where the effects of the pin can 
be treated perturbatively is simplest to analyze, as was 
first done in Ref. 35. This can be done in close analogy 
to Sec. IV B, extending it to a finite h. The governing 
Hamiltonian is given by 

n = dxdz [K{d^uf + B{d^uf] 

~v j dzcos[G{u{0,z) + hz)], (4.31) 

with u{x,z) fiuctuating around and to zeroth order 
in the pin strength v exhibiting translationally invariant 
correlations: 



{{u{x,z)-u{0,0)f)o 



■ks/KB 



In 



B 



(4.32) 



F,i(r) « E,i{T)-kBT In 



C(T)' 



To the same zeroth order the average vortex density is 
uniform, uq. To compute it to the lowest nontrivial order 
in w, we expand the Boltzmann weight to first order in the 
pinning potential, v and utilize the phonon correlations 
of u{x, z), Eq.4.32 about the tilted state. This leads to: 
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ncv . 



{n{x, z)) — no ~ - — —Re 



Rc 



,,i.Gx 



dz' 



r,iGh{z—z') 



^^gr,{\x\/a)e-\^\/'^cosi2nnox), 



(4.33) 



with the transverse field h giving an exponential fall off of 
Friedel oscillations on scale A^, Eq.3.51,'^^ and preexpo- 
nential function gr,{\x\/a) approximated by a power-law 



5r,(k|/a) 



:,(^) , a«|x|«A^, 

(4.34) 



c^, dri 0(1) dimonsionlcss constants. 

The exponential decay of density correlations, Eq.4.33 
vindicates a perturbative treatment (in v) at large trans- 
verse field and temperature, where the decay length 
is the shortest scale in the problem. 

(b) Commensurate (pinned) transverse Meissner state, 
h < hci{T): 

Vortex correlations are also simple to analyze in 
the low tilt-field, low-temperature commensurate h < 
hci state, where, because of the commensurate- 
incommensurate (tilting/roughening) transition we ex- 
pect a behavior that is qualitatively different from that of 
the incommensurate state. We expand the Hamiltonian, 
Eq.4.31 in vortex lattice displacements w{x, z) about the 



J 



nontrivial zero-temperature distortion uq{x,z), Eq.3.48 
(with Us = 0) characterizing this state: 



u{x, z) = uq{x, z) -h w{x, z). 
To quadratic order in w we find 



(4.35) 



(5H = / dxdz 



|(a.«;)^ + |(5.«;)^' 



1 



+ -V (27rno) / dzw(0,z)^. 



(4.36) 



where we used the fact that in the commensurate state 
wo(0, z) = 0. Clearly then, because the pinning nonlin- 
earity is localized at a; = 0, vortex correlations of w{x, z) 
on top of the ground state background uo(a;, z) are iden- 
tical to those of the /i = state. Namely, on scales longer 
than the pimiing length ^, w{x, z) fluctuations are that 
of a field that is harmonic in the bulk, but pinned at the 
X = Q boundary. This identification allows us to take 
over results from Sec. IV C. Using the simplified version 
for uo{x, z), Eq.3.49, we therefore find the average vortex 
density 



{n{x,z)) = no - nodxUoix,z) +2nG^i.e 



no - no\/ -^Y^hz sgn(a;)e~v btNI + 2nG ( ^73 j cos [27rno(a; -|- uo(a;, ^;))] , 



(4.37) 
(4.38) 



that exhibits power-law decaying Friedel oscillations 
(non-perturbative in v) as well as a smooth compression 
and dilation of the vortex lattice around the pin for z ^Q. 
Both features can be clearly seen in the vortex configu- 
ration displayed in Fig. 13. 

(c) Incommensurate (soliton-tilted) state, hci(T) < 
h^hc2{T): 

We now turn to the intermediate incommensurate 
regime of h and T, with h > hciiT). Because there 



is only a single CI transition (that across hd, with hc2 
simply a crossover), we do not expect any qualitative 
change in correlations from those found above for large h 
and T . However, examining the expression for {n{x, z)), 
Eqs.4.33, 4.34, the prefactor (a/A^)''/2"\ that is large 
for ?7 < 2 and A^ » a, suggests a failure of the per- 
turbative treatment for T < Tp and for a sufficiently 
low h. Indeed, as illustrated in Fig. 17, in suc h a regime 
A^ will exceed the pinning length iy/BjK. Since for 
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T < Tp, on scales longer than BjK, even a (bare) 
weak pinning potential becomes comparable to the elas- 
tic energy and therefo re can not be treated perturbatively. 
Consequently, for S^^/bJK < there opens up an in- 
termediate regime, S^^B/K < x < A^, where above 
perturbative (in v) analysis leading to exponential decay 
of Friedel oscillations, Eqs.4.33, 4.34 fails. The crossover 
boundary (T) separating the perturbative and nonper- 
turbative regimes is clearly given by f^\J BjK = A^, 
which (not surprisingly, but reassuringly) is equivalent 
to K{T) « /ic2(T) « a/27rC(r) found in Eq.3.36. 




nonperturbative 



FIG. 17: A schematic of the length-scale dependent pin- 
ning coupling v{x), indicating its intermediate nonperturbar 
tive regime for h < hc2{T). 

As schematically illustrated in Fig. 17, in the parlance 
of RG, for T < Tp the dimensionless pinning coupling 
v{x) at short scale x grows as x^~^^^. For h > hc2{T), he- 
fore v{x) reaches a nonperturbative regime (i.e., reaches 
0(1)), this growth is "quenched" by the tilt field h on 
scale longer than A^, beyond which v{x) ~ e""^/^^ de- 
creases exponentially, validating the fully perturbative 
treatment of case (a) above. In contrast, for h < hc2{T), 
v{x) becomes large on scale £,^B/K < A^, requiring a 
nonperturbative analysis on longer scales. 

Above discussion clearly indicates a need for separate 
treatments of three distinct regimes of length scales (i) 
a <x < £,^/B/K, (ii) i^jBjK < x < Xh, (iii) A^ < x. 



In the regime (i) of shortest lengths scales, the pinning 
potential in 7Y, Eq.4.31 can be treated perturbatively, 
with average density given by the x <C A^ limit of the 
result in Eq.4.33, with exponential factor approximately 
1. On longer scales, in regime (ii) the pinning potential 
is comparable to the elastic energy and vortex configura- 
tions displays a soliton array of density ns{h), illustrated 
in Fig. 15. In this regime ^y^B/K < x < X'^ < Xh seg- 
ments of vortex lines on the corresponding scales shorter 
than appear to be strongly pinned in the vicinity 
of the columnar defect. Hence on these scales we ex- 
pect vortex correlations of the transverse Mcissner state 
computed in Eq.4.38, with only small corrections from 
soliton array fluctuations. Finally, on scales z > nj^ 
(corresponding to x > X^ of regime (iii)), it is clear from 
Fig. 15, that vortex lattice displays an average tilt ans{h), 



u{x, z) = angZ + w{x, z), 



(4.39) 



with small fluctuations w{x,z) about the soliton state 
arising from vibration of the soliton array. We expect 
the corresponding effective Hamiltonian for w{x, z) to be 
given by 

Sn = dxdz {K{d^wf + B{d^wf] 

~v / dz cos {G{w{0, z) + augz)), (4.40) 



with K K. K, B K. B the effective moduli of the soliton- 
tilted vortex lattice and v w v{S,/ a^)^^/"^ a weak pinning 
potential, reduced by thermal fluctuations from scales 
a2 < ^; < ^ of regime (ii). On longer scales z > 
(regime (iii)), the soliton spacing is the shortest scale 
in the problem (beyond the lattice spacing), and analysis 
of w correlations using 5TL can be done perturbatively in 
V similar to that of regime (a), that led to Eq.4.33. As 
can be seen from comparing Hamiltonians in Eq.4.31 and 
Eq.4.40, the main qualitative difference is the replace- 
ment of h by ans{h) inside the pinning potential, that 
leads to the corresponding replacement of A^(/i) by the 
longer length Xh{h) in the exponential decay of Friedel 
oscillations of the average vortex density. Putting these 
results together for /id < h < hc2 we predict 



{n{x, z)) — no ~ cos(27rnoa;) < 



npyg j K 



\x\ 

- \ v/2 



'7-1 



a < |x| < i^B/K, 



B \\x\) 



ngvXh 



77-1 



-\x\/\h. 



i^B/K <\x\<Xh, 
Xh < \x\. 



(4.41) 



V. FINITE DENSITY OF INDEPENDENT 
COLUMNAR DEFECTS 



tion, no genuine tilting CI phase transition is possible 



So far we have focused on an idealized problem of a 
single columnar defect. As discussed in the Introduc- 
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in this case, since pinning energy density and the asso- 
ciated lower critical field hd, Eq.3.22 vanish in the 2D 
bulk thermodynamic limit. In this section we extend our 
results to a physically more interesting case of a finite 
dilute concentration 1/d of columnar defects. 

A full treatment of such a highly nontrivial problem 
is beyond the scope of the present paper and has been 
a subject of numerous studies. ^^'^^'^^"^^ Here we will be 
content with an approximate analysis of a dilute con- 
centration of columnar defects on intermediate length 
scales, where vortex lattice response around each pin can 
be treated independently. However, wc expect that for a 
fixed defect concentration (even if dilute), the system will 
crossover to the anisotropic (Bose) vortex glass collective- 
pinning phenomenology and our results will break down 
on sufficiently long scales. This is in the spirit of other 
problems in physics, most notably the Kondo effect^'', 
where for a sufficiently dilute concentration of impuri- 
ties, on intermediate scales local moments can be treated 
independently, but may order magnetically at sufficiently 
long scales and low temperatures. 

The main qualitative effect of a finite pin spacing d 
should be clear from the analysis of Sec. HID, in partic- 
ular from Eqs.2.10, 3.48 and vortex configuration illus- 
trated in Fig. 15. There wc have shown that a vortex 
lattice distortion of waveleng th Az along the defect pen- 
etrates to length A_L = Xz\/ B/K /2'k along x into the 
bulk, away from the defect. For a finite concentration 
of columnar defects, we therefore expect long wavelength 
elastic distortions with > K/ B, corresponding 

to overlapping distortion clouds of neighboring pins to be 
cut off at scale d along x. 

To demonstrate this in detail, we need to generalize our 
results in Eqs.2.10, 3.48 to a nontrivial boundary condi- 
tion on the clastic distortion u{x, z) at x = d, and use it 
to recalculate predictions of previous sections, most im- 
portantly the energies of the commensurate (aligned) and 
incommensurate (soliton) states. The correct qualitative 
physics that we are after here can be obtained by the use 
of periodic boundary condition on u{x, z) with period d, 
or even more simply (and a bit cruder) a boundary con- 
dition of simply cutting off x integrals beyond length d/2 
around each pin, corresponding to the Dirichlet bound- 
ary condition on u{x, z). 

To see this in more detail, wc recalculate the energy 
per columnar defect E = 'H[uo{x, z)], Eq.2.5, 

1 fL /•d/2 

E = - dz dx [K{d,uo - hf + B{da,UQf] 

^ Jo J -d/2 

r-L 

dz COS {Guq{Q,z)), (5.1) 

with the above boundary condition at \x\ = d/2 and 
uo{x, z) given by Eq.3.48. Since Euler-Lagrange equation 
ensures that contributions from each of the three terms 
balance each other, E approximately reduces to 

l.d/2 



(5.2) 



Calculating above expression wc find, that, for L < 
2-Kd\jKlB it reduces to the previously calculated single- 
pin result given in Eq.3.32. In the opposite regime, 
L > 2'Kd^KIB of interest to us, E crosses over to 



E{d) « Ec{d) ^Kad{hd-K)]<S, + \Ka^yNl, 



with bulk energy per pin 

Ec{d) = Kh^Ld-vL, 

leading to 



and 



hci{d) 




KAnd 

as quoted in the Introduction. 



In 




(5.3) 

(5.4) 

(5.5) 
(5.6) 



E 



Jo Jo 



dxK 



dzUo{x,z) - hese 



vL. 



VI. CONCLUSIONS 

In this paper we studied a finite-temperature re- 
sponse of a planar vortex array to an in-plane tilting 
of an external magnetic field away from a dilute con- 
centration of pinning columnar defects. We found that 
the vortex lattice tilting proceeds via an interesting fi- 
nite temperature, finite transverse field commensurate- 
incommensurate transition at /ici(T), driven by a prolif- 
eration of solitons, as illustrated in the phase diagram. 
Fig. 3. A sensitive dependence of this lower-critical field 
(that vanishes for a single pin) on the columnar defect 
spacing should be experimentally testable by varying the 
heavy-ion irradiation flux used to create pinning tracks. 
We show that at low-temperatures, for h > hci{T) the 
vortex array exhibits a highly nontrivial soliton-like dis- 
tortion as a compromise between inter- vortex interaction, 
and the pinning and diamagnetic energies. We show 
that this nonlinear response persists up to a large upper- 
critical tilting angle, tan^c2 ~ a/^{T), beyond which sys- 
tem recovers a full linear transverse susceptibility with 
B± « H±. We expect that these and many other de- 
tailed predictions should be directly testable in computer 
simulations. Although much more difficult, it is our hope 
that the theory presented here can be furthermore tested 
in mesoscopic samples of artificially layered superconduc- 
tors. 

More importantly, we expect that a number of these 
features will carry over to a transverse field response in 
bulk superconductors. Extending our two-dimensional 
planar results and exploring their impact on phenomenol- 
ogy of bulk samples remains an important and challeng- 
ing problem. 
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APPENDIX A: BOUNDARY HAMILTONIAN VIA 
A FUNCTIONAL INTEGRAL 

Problems where it is possible and convenient to elimi- 
nate bulk degrees of freedom, thereby reducing the prob- 
lem to a lowcr-dimcnsional one arc quite common in 
physics. In this appendix, for completeness, we present a 
general functional-integrals approach to such problems, 
and apply it to the problem of a planar vortex lattice 
pinned by a single columnar defect treated my more 
pedestrian methods in the main text. 

To this end, we consider a field (/'(x, z) defined on a 
coordinate space r = (x, z) that we split into bulk x and 
boundary z subspaces. The energetics is governed by a 
Hamiltonian and the partition function is given by 
a standard functional integral over the field ^{r) {ksT = 
1): 



-n[<i>] 



(Al) 



Z can be cquivalently expressed as an integral over field 
(/)(r) constrained on the boundary x = to be (j)(0,z) = 
(/>o(z), followed by an integral over the boundary fields 
(/>o(z). Explicitly, the former is implemented by a func- 
tional i5-function, i5[<?f>(z)] = riz^C"!^!^)) giving 

Z = J [#o(z)#(x,z)] 5[</)(0,z) - 0o(z)] e-^[*l, (A2) 

= /[#o(z)] e-^°l*°(")l. (A3) 



The effective boundary Hamiltonian Ho[4>o{z)] defined 
above can be expressed using the Fourier representation 
of the functional 5- function, leading to 



J [dcj>dX] e-«[^l+'/.M-)(<^(o,z)-^o(z))_ 

(A4) 



For a general Hamiltonian 7i[0], above computation can 
only be performed via a formal cumulant expansion. 
However, for a special case of a quadratic Hamiltonian 



defined by a correlation function G{r) (with G~^{r) its 
inverse), all cumulants reduce to a power of G(r), i.e., 
obey Wick's theorem, equivalent to the Gaussian integral 
identity 



f 

J — ( 



^+\4> - 



2tt\ 
a J 



1/2 
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(A6) 



A Gaussian functional integrations over (f>{r) and A(z) 
(dropping an inconsequential (/)o-iudependent constant) 
then gives 



-i ^, A(z)G(0,z-z')A(z')-i A(z)<^o(z) 



-i/,,, 0o(z)G-i(O,z-z')0o(z') 



(A7) 



Applying this simple result to the (l-l-l)-dimensional 
(planar) vortex lattice pinned by a columnar defect at 
X = leads to 

no[uo{z)] = i / uo{z)G^\z - z')uo{z'), (A8) 



with 



(A9) 



an inverse of the bulk propagator G(0, z) evaluated at 
X = 0. In g^-Fourier space the latter is easily evaluated 



G{x = 0,q,) 



dqx 



1 



27r Kcg + Bql' 
1 1 



and leads to 



G^Hq,) = 2VBK\q, 



(AlO) 
(All) 

(A12) 



thereby confirming the resulting for Hoiuoiqz}], Eq.2.12 
obtained in Sec. II B by a different method. 



APPENDIX B: BULK CORRELATION 
FUNCTIONS IN A "PINNED" STATE 

In this appendix, as a model of phonon correlations 
in the "pinned" vortex state, we study bulk correlation 
function of a Gaussian field constrained at r = (0,z) 
by a "massive" boundary Hamiltonian 'Hp[(f>{x,0)]. As 
argued in the main text, on sufflciently long scales the 
boundary Hamiltonian can simply be implemented as a 
hard constraint on the field to vanish at r = (0,z). To 
this end we compute the asymptotic generating function 

Z[j{r)] = y[d0(x,z)]J[</.(O,z)]e-^[*l+X^'W^«, (Bl) 

from which, by differentiation with respect to j(r) all 
n-point correlation functions of (p{r) can be obtained. 
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As in Appendix A, we have implemented the Dirichlet The kernel r[q, q'] is given by 
boundary condition on 0(r) via a functional 5-fimction. 
Representing the latter in its (functional) Fourier form 
as in Eq.A4, using harmonic Hamiltonian, Eq.A5, and 
performing a Gaussian integral over ^(r), we find 

Z\j{r)] = eiL.>=^^^^M^^(^'\ (B2) 



r[q,q'] = (27r)'^^(5'^^(q, +q'JG(q,,q,) 



(27r)'^-J''-(q.+q;)-G(q;,-q.) 



1x 



G(q;',q.) 



(B3) 



and leads to real-space 0(r) 2-point correlation function 

Gpi„ned[x,x';z] = ((/)(x,z)(/)(x',0))pin„ed, (B4) 

that in q^-Fourier space is given by 

G(x,q,)G(x',-q,) 



Gpinned[x,x';q^] = G(x-x';q^) 



G(0,qz) 



(B5) 



In the special (l+l)-dimensional case of G ^{qx, Qz) 
Kq^ + Bq^ of interest to us in the main text, 



G{x,qz) = J 



dqx 



e-(f)'l9»lkl 
2y/BK\qz\ 



(B6) and leads to 



(B7) 



Gr,\ 



pinned {-^t ^ ; 



1 



2VBK\qz^ 
G{\x-x'\,z)-G{\x\ + \x'\,z). 

I 



(B8) 
(B9) 



utilized in Eq.4.23 of the main text. 

A simpler way to derive above "pinned" correlation 
function result is to note that the boundary condition 
0(0, z) = is automatically explicitly satisfied by the 
odd part of (f){x,z). It is also satisfied by a subset of 
(an independent) even part of (/)(x, z) that vanishes at 
X = 0. Naively, such even/odd field decomposition does 
not represent a simplification since a constraint on the 



even part must still be enforced. However, it is clear that 
for correlations on the same side of the pin, a constrained 
even part of the field has identical correlations to that 
of the odd part of the field. Hence, Gpinnod[x, x'; z] for 
xx' > (i.e., on the same side of the pin) is simply given 
by twice the correlator of the unconstrained odd part, 
(/)(x,z) 5(0(x,z) - 0(-x,z)) 



J 



Gpinned[x,x';z] = ^ ((</.(x, z) - </.(-x, z)) (,^(x' , 0) - ,^(-x' , 0)))o, (BIO) 
= (</.(x,z)</.(x',0))o- Wx,z)0(-x,O))o, (Bll) 

I 



giving the result in Eq.B9. 



We conclude this appendix with a computation of the 
generating function Z[j(r)], that extends above analysis 



22 



to scales shorter than the pinning length. To this end, The corresponding generating function is then given by 
we supplement the harmonic bulk Hamiltonian with a 
pinning one at the boundary, 7ip[(p{0, z)] with the full H 
given by 

= \j </.(r)G-i(r-r')</.(r')+ j ^''»(x)Wp[,/.(0, z)], 

(B12) 



ZW)] = j [#o(z)]e-«''[*°Wl j [#(x,z)] 5[4>{0,2,) - </.o(z)] e 



0(r)G-i(r-r')0(r')+ /_.j(r)0(r) 



(B13) 



r 



As in Appendix A, representing the functional (5-function with 
in its Fourier form, integrating over the harmonic bulk 
and boundary Fourier fields ^(r), A(z), we find 



ZW)] 



[d(/)o(z)]e' 



-W[<^o(z),j(r)] 



(B14) 



W[<^o,i] = Hp[U^)]-\j^ j{v)G{v-v')j{v') 

+ ^ [ M^)- I Go(xi,z - zi)j(xi,zi) 

^ J z,z' L J Zi ,Xi 



(B15) 

G-1(0,Z-Z') M^')- f Go(x2,z'-Z2)i(x2,Z2) 

I 



For j{r) and HpiM^)] = J^cosM^), W[</'o(z),0] 
simply reduces to the Hamiltonian, Eq.2.11 localized on 
the pin, with the integrated out bulk degrees of free- 
dom reflected in its long-range elasticity. On the other 
hand, for a finite j(r), but strong pinning, on sufSciently 
long scales e^''^p['*"(^)l simply acts as a hard constraint 
(/)o(z) = 0, reducing 2'[j(r)] to the previously found re- 
sult given in Eqs.B2, B3, B5. 

In the pinned (commensurate) phase, a columnar de- 
fect pins a single vortex line, corresponding to a "con- 
finement" of field (poiz) to a single minimum of the 



cosine and allowing us to approximate 7ip[(f>o{z)] = 
—V /^cos0o(z) ~ const. -I- |w /^(^o(z)^ by a harmonic 
"spring". We can therefore integrate over ^o(z) in 
Eq.B14, obtaining (up to an unimportant multiplicative 
constant): 



Z\j{v)] 



gi jWrpin„ed[r,r']j(r')_ ^^^^^ 



with 



J 



rpi„„ed[r,r'] = / G(x,z;-z)G-i(0,zi-z;) [t;<5(zi-Z2) + G-i(0,zi-Z2)] ' G(x',z'2 - z')G-i(0,Z2 - z^) 
+ G(r-r')-/" G(x,z-zi)G-i(0,zi-Z2)G(x',z'-Z2). (B17) 

>^Zi,Z2 



In Fourier space, this becomes 
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rpinned[q,q'] = (27r)'^^ 5"^^ (q, + q'J [(27r)'^« ^-^^ (q, + q'JG(q, , q,) 
- G(q,,q,)G-i(x = 0,-q,)G(q',,-q 



r + r;-i(x = ().q,)- 



(B18) 



For a system that is translationally invariant along z, 
G^^(x = 0,qz) generically vanishes at long wavelengths 
q^ ^ 0, and the result reduces to that of a hard con- 
straint, given in Eqs.B2, B3. In more detail for the 
(l-l-l)D vortex problem at hand, G~^(x = 0, q^) = 
{a/2TT)'^2y/KB\qz\, showing that, as asserted in the main 
text, the crossover to the hard constraint happens on 
scales longer than the pinning length, ^ ^. 



be computed as (minus) the imaginary part of Hilbert 
transform of 



with the latter function ana- 



lytic in the upper-half plane. As a side benefit the real 
part of H[^] gives 



Hf 



y2 + i 



(C8) 



APPENDIX C: HILBERT TRANSFORM BASICS 



APPENDIX D: SOLITON SOLUTION OF THE 
SINE-HILBERT EQUATION 



In this appendix, for completeness we summarize some 
of the basics of Hilbert transforms necessary to derive 
results in the main text and in the Appendixes. 

Hilbert transform (f){y) = Il[(p{x)] of a function 4>{x) is 
defined by 



m(x)] = -P / dx- 

TT J_oo X 



y 



(Cl) 



where P stands for the principal value of the integral i.e., 
with the singular point x ~ y excluded. 

Hilbert transforms of standard functions can be usually 
computed by relating it to a contour integral in a complex 
plane. For example, Hilbert transform of sinx and cos.x 
can be computed as real and imaginary parts of Hilbert 
transform of e'^ 



1 r 

H[e"] = -P / 



dx- 



x-y 



(C2) 
(C3) 



with last expression obtained easily by contour integra- 
tion, taking advantage of analyticity of in the upper- 
half plane. Above result then leads to 



H[sina;] = cosy, 
H[cosa;] = — siny. 



(C4) 
(C5) 



More importantly for the problem of the vortex lat- 
tice at hand, we compute the Hilbert transform of a 
Lorentzian 



Hf 



'x^ + V 



1 



_P / dx 2X7' 

TT 7_oo x-yx^ + 1 

-y 



2/2 + 1' 



(C6) 
(C7) 



by noting that it is related a semi-circular contour inte- 
gral over C in the upper-half plane. Equivalently, it can 



In this appendix we verify that the soliton solution 

1 



0s (z) = -2ArcTan- (Dl) 
z 

indeed satisfies the Euler-Langrange integral equation 

,(t>s{z) - Mz') 



dz' 



{z - z'f 



+ sin 0s (2) = 0, 



(D2) 



where from now on, all the integrals are understood in 
the sense of a principal part, a physically dictated regu- 
larization. To this end, using a relation P / ^^^^,^5 = 
and integrating by parts, we note that the sine-Hilbert 
equation can be rewritten as 



Now using 



-9.H[0s(z')]+sin0s(z) = 0, 
-H[5^-0s(/)]+sin0s(z) = 0. 



(D3) 
(D4) 



(D5) 



and Hilbert transform relation from Appendix C, Eq.C7, 
we find 



-2z 



z'^ + 1 



Then calculating 



sin (ps{z) = —sin 
-2z 



2ArcTan- 

z 



(D6) 

(D7) 

(D8) 



shows that indeed 0s (2;) satisfies the sine-Hilbert equa- 
tion, Eq.D4. 
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APPENDIX E: SINGLE SOLITON ENERGY 

In this appendix we compute the energy Ei[h] 
'Ho[(f)s{z)] of a single sohton, 



-2ArcTan- 



1 



Z- Zo 

for the sine-Hilbert model, defined by a Hamiltonian 

/to — /to + /to, 



(El) 



= — I I dzdz' 
An , 



(z) - 4>{z') - h{z - z') 



z — z' 



-I 



dz COS (p{z), 



(E2) 



where all integrals are implicitly understood to range over 
the system size, with —L/2 < z < L/2. Expanding the 
square of the elastic part, Ef = 'HQ[(j)s{z)], we find (to 
accuracy of 0{1) for L oo): 



Ef = 



1 

2^ 



dzdz' 



(psiz)"^ - (l>s{z)4>s{z') _ ^ Mz) - 4>s{z') 



^ j dzcj>,{z)B.[d,(f>s{z)]-h J dzR[Mz)] + ^h^L\ 



2 f 2 



4tt 



-h^L 



(E3) 
(E4) 



These integrals can be computed utilizing Hilbert trans- 
forms worked out in Appendix C. Using Eq.C7, the first 
term, Ef^^ can be integrated by parts 



= 27rln- 



dz- 



P + 1 



-ArcTan- 



(E5) 
(E6) 

(E7) 



Less formally, this elastic contribution can be computed 
by going back to the expression 

2 



dzdz' 



(ps{z) - <i)s{z') 



(E8) 



and noting that because (t>s{z) vanishes for z ^ 1 and 
equals 27r for i < — 1, finite contributions to the elastic 
energy arise only from regions (^ < — 1, 5' > 1) and {z > 
1, z' < —1). To accuracy of 0{1), this reduces the soliton 
elastic energy to 



dz 



L/2 



l\2- 



(E9) 
(ElO) 



in agreement with the more formal analysis above. 

The second contribution to Ef in Eq.E4 can also 
be computed by integrating by parts and noting 
that H[^j,(£)] (that can be explicitly computed giving 

H[(/)s(z)] = — In ) vanishes at the boundaries 

\_(L/2) +1J 

of the system, z = ±L/2 



Ef, = -h 



L/2 



L/2 



dzHlMz)], 



(Ell) 



= h dzm[di(l)s{z)], (E12) 

J-L/2 

= -/i(2L-27r). (E13) 

A single soliton pinning contribution Ef = H^lcpsiz)] 
is also straightforward to calculate using solution (ps{z), 
Eq.El. We find 

Ef = - dz cos (j)s{z), (E14) 

J-L/2 
rL/2 

= -/ d£(l-2sin2(</)^/2)), (E15) 

J-L/2 



.L/2 




/ . dz 




/-L/2 





= -L-\-2-K. 



(E16) 
(E17) 



Combining above contributions inside Eq.E2, we ob- 
tain the expression for a single soliton dimensionless en- 
ergy 



1 



eL 



El = —h^L^ -1 + 2^171— -h{2L-27T), (E18) 



47r 

—h'^P - L + 2-K\nL-2hL 
4n 



(E19) 



used in the main text. 



APPENDIX F: TWO-SOLITON ENERGY: 
INTERACTION 

In this appendix we give a few technical details for 

the computation of the soliton interaction energy. For 
two far-separated solitons (of interest in a dilute soliton 
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approximation, valid for h <C hc2) with |ii — Z2I 3> 1, wc 
can approximate the exact two-sohton solution by a sum 
of two one-soliton solutions 

<l>2s{Zl,Z2) ~ (paizi) + (ps{z2)- (Fl) 



The two-soliton interaction is then determined by two- 
soliton energy 



E2{Zl,Z2) « Hol^sl + (I>s2], (F2) 

1 f ( (ksl{z)+ct>s2{i)-<f>sl{z')~<l>s2{z')-Hz-z') \' /•,, , 

^ An \ g _ I - / dz cos[(f)si{z) + (I)s2{z)], (F3) 

= Ec + 2Esi+Vs{zi- Z2), (F4) 



I 

that consists of the zero-soliton contribution Ec, Eq.3.20, Vs{zi — Z2) = V^{zi — Z2) + VP{z\ — Z2) given by 
two one-soliton contributions Esi-, and soliton interaction 



dzdz 



,(j)sl{z)(j)s2{z) - (I>sl{z)(ps2{z') 



(z - z'Y 

Vg^{zi-Z2) = / [1 - COS ^si(z) 008(6,2(2) + sin (6.,! (z) sin c/),,2(z)] 



(F5) 

(F6) 



Manipulations similar to those for the computation of tion gives 
a single soliton energy give 



^{zi-Z2) ~ ~y dZ(t)sl{z)B.[dB(l)s2{z 

Z — Z2 



(F7)K''(^i-^2) = 2 



dz 



(2z - zi - 



((Z-Zi)2 + I)((z-Z2)2 + 1) 



, (Fll) 



4 / dz 



(Z - Z2f + 1 



ArcTan- 



Z — Z\ 



-, (F8) 



47rln- - 2 
2 



47r In - - 27r In \(zi - Z2f + 4 
2 



In \(z - Z2f + 1] , . 
2 . ^^ (Fio) 



= 47r, (F12) 
that, together with Eq.FlO gives the soliton interaction 

\Ll2f 



27rln 



z2 + 4 



+ 47r, 



(F13) 



with a simplifying approximation above valid for |zi,2| <C 

i/2 00. Similarly, a straightforward contour Integra- used in the main text. 
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